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We investigate non-classical effects such as fractional revivals, squeezing and higher-order squeez- 
ing of photon-added coherent states propagating through a Kerr-like medium. The Wigner functions 
corresponding to these states at the instants of fractional revivals are obtained, and the extent of 
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I. INTRODUCTION 

Non-classical effects such as revivals and squeezing 
exhibited by quantum wave packets while propagating 
through different media are of great interest in the con- 
text of quantum information processing with continuous 
variables[l|. Both the nature of the medium in which 
the state propagates and the precise initial state consid- 
ered play a crucial role in determining the subsequent 
dynamics of the state. With the experimental realiza- 
tion of several non-classical states in quantum optics 0, 
it has become important to examine more closely their 
dynamical properties. Detailed investigations have been 
carried out, for instance, on the revival phenomena dis- 
played by an initial coherent state evolving under dif- 
ferent HamiltoniansQ- In particular, the non-classical 
features displayed by a wave packet which is initially a 
coherent state while it propagates through a nonlinear, 
Kerr-like medium have been examined. Its squeezing and 
higher-order squeezing properties have been discussed Q, 
and signatures of wave packet revivals and fractional re- 
vivals have been shown to be captured in the expectation 
values of appropriate observables Q . 

In this paper, we examine the effects of the departure 
from coherence of the initial state on the subsequent non- 
classical features it displays while evolving in a Kerr-like 
medium. In order to be able to compare with exacti- 
tude these effects with those that arise in the dynamics 
of an initial coherent state (CS), it is essential to con- 
sider such initial states as display a precisely quantifiable 
departure from coherence. The photon-added coherent 
states (PACS) are promising candidates for this purpose 
@. We examine a spectrum of non-classical phenom- 
ena such as revivals, fractional revivals, squeezing and 
higher-order squeezing displayed by an initial PACS as 
it propagates in a Kerr-like medium, and identify the ef- 
fects arising from the departure from perfect coherence 
of the initial state. 

The plan of the paper is as follows: In the next sec- 
tion, we review briefly the salient features of revivals and 
fractional revivals of a wave packet, and identify quan- 



tifiers which provide signatures to distinguish between 
fractional revivals of initial wave packets that are ideal 
CS from those that are PACS. In Section III, we exam- 
ine the dependence of squeezing and higher-order squeez- 
ing effects on the extent of coherence of the initial wave 
packet. In the final section, we estimate the degree of 
non-classicality displayed by a CS and by a PACS at 
specific instants while evolving in a Kerr-like medium, by 
quantifying the extent of negativity of the corresponding 
Wigner functions. 



II. WAVE PACKET REVIVALS 

A wave packet propagating through a nonlinear 
medium may display revivals and fractional revivals at 
specific instants of time. In the one-dimensional case, if 
the initial wave packet 1^(0)) lS a superposition of a suf- 
ficient number of basis states, it is expected to revive at 
periodic intervals of time^though these intervals could 
be quite long in practice 7]. A full revival at time T 
is signaled by the fact that the autocorrelation function 
\(ip(0)\ip(T)}\ 2 returns to its initial value of unity. 

In a generic nonlinear medium, for wave packets which 
are peaked sufficiently sharply in energy about some 
eigenvalue labeled by uq, we may expand the energy spec- 
trum E n in a Taylor series about E no , and retain only 
terms up to the second order in (n — no) . It can then 
be shownQ that fractional revivals of the wave packet 
occur at certain instants of time between successive re- 
vivals, when the wave packet splits into a superposition of 
a finite number of spatially separated sub-packets, each 
of which closely resembles the original one. This non- 
classical phenomenon arises due to very specific quantum 
interference properties between the basis states compris- 
ing the original wave packet [To| . The dynamics of a wave 
packet evolving in a Kerr-like medium is governed by the 
model Hamiltonian 



H = h\a 



where N = a^a and a , a 
creation operators, so that 



= h X N(N-l), (1) 
' are photon annihilation and 



[OifflT] = 1. The positive 
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constant \ represents the susceptibility of the nonlinear 
medium. (The same Hamiltonian is applicable to a Bose- 
Einstein condensate propagating in a three-dimensional 
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optical lattice. In this case, a and a) are atom anni- 
hilation and creation operators, and x characterizes the 
energy needed to overcome the inter- atomic repulsion.) 

For ready reference, we first state the results^ for an 
initial coherent state \ip(0)) = \ot), where a\a) — a\a) 
and the complex number a = (x +ipo) / \[2 labels the CS. 
Such an initial state can be shown to revive periodically, 
with a revival time 

T rev = 7t/ X . (2) 

Between t = and t = T rcv , fractional revivals oc- 
cur at times t — Trj/(kx), where k = 2,3,... and 
j = 1, 2, . . . , (fc — 1) for a given value of k. At time 
(j/k)T Tev , the initial Gaussian wave packet splits into k 
spatially distributed sub-packets that are similar to it- 
self. Full revivals of the wave packet can be discerned 
by observing the time evolution of the expectation val- 
ues (first moments) of the operators x = (a + a))/y/2 
and p = (a — al)/(iV2), while distinctive signatures of 
the fc-sub-packet fractional revivals manifest themselves 
in the temporal behavior of the k th moments of x and p. 
The time dependence of all moments of x and p can be 
obtained from the general result 



= a s v r e^ u ^^ cos 2sxt ^ 

3 — ix[s(s— l)+2rs] t—iv sin 2s\t 



x e 



(3) 



where r and s are non-negative integers, and v = |a| 
is the mean number of photons in the CS. Writing 
(ip(t)\ x \ip(t)) — (x(t)), etc., it is evident that (x(0)) = x n 
and (p(0)) = po . We find 



<*(*)> 



-v (1 — cos 2\t) 



\xq cos {y sin 2\t) 



p sin (z/ sin 2x*)]. 



(4) 



with a similar expression for (p(t)). Owing to the ex- 
ponential factor on the RHS we find that, for sufficiently 
large values of v, (x(t)) remains essentially equal to (x(0)) 
for much of the time, except for significant and rapid vari- 
ation near instants of full revival. This behavior is shared 
by (p(t)). A sudden change in (x(t)) or (p(t)} thus sig- 
nals the occurrence of a revival of the initial wave packet. 
Similarly, writing (ip(t)\ x 2 \ip(t)) = (x 2 (i)), we can show 
that 



2(x 2 (t)) 



1 



-v (1 — cos 4xt) 



x [Oo - Po) cos ( 2 Xt + v sin A X t) 

+ 2x p sin (2\t + v sin 4%£)] . (5) 



Once again, for large values of i>, (x 2 {t)) is more or less 
static around the value (x 2 (0)) except at times close to 
7r/(2x) = |T rov , when the 2-sub-packet fractional revival 
occurs. Similar behavior is displayed by (p 2 {t)). An ex- 
plicit calculation reveals that, in general, the signature 
of a fc-sub-packet fractional revival is mirrored in the dy- 
namics of (x fc (i)) (equivalently, of (£))). It follows 



that the plot of Ax (the standard deviation in x) versus 
Ap (the standard deviation in p) captures the occurrence 
of a 2-sub-packet fractional revival. Likewise, the square 
of the skewness in x (or p) mirrors the occurrence of a 
3-sub-packet fractional revival, while the temporal be- 
haviour of the kurtosis of x (or p) enables us to detect 
the appearance of a 4-sub-packet fractional revival. 

Next, we turn to the case in which the initial state 
is a PACS rather than a (perfectly coherent) CS. The 
normalized m-photon-added coherent state \a,m) (m = 
1,2,.. .) is defined as 



\a, m) 



(at) m | a) 



_ (at) m |q) 

y/(a\a m at™ \a) ^m\L m {-v\ 



(6) 



where L m is the Laguerre polynomial of order m. The 
PACS just defined is also a nonlinear coherent state in the 
sense that it is an eigenstate of a nonlinear annihilation 
operator, namely, 



1 



1 + at a 



a \a, m) = a \a, m) 



(7) 



A PACS of this sort can be produced in laser-atom in- 
teractions under appropriate conditions 6] . Clearly, with 
increasing m (the number of photons 'added' to the ideal 
CS | a)), there is increasing departure from perfect co- 
herence. It is straightforward to see that, when an initial 
state \a,m) evolves under the Hamiltonian in Eq. 
revivals and fractional revivals of the wave packet occur 
at the same instants as in the case of a coherent initial 
state, and that their appearance is mirrored in the time 
evolution of the expectation values of the different mo- 
ments of the quadratures x and p, as before. However, 
the actual dynamics of these expectation values is much 
more complex, even for small values of m: even a small 
departure from Poisson number statistics in the initial 
state leads to significant changes in the evolution of the 
wave packet. 

To illustrate this explicitly, we find first the counter- 
part of Eq. for the general initial state \ip(0)) = 
la, m). The result can be written in closed form as 



(a tr a r+s 



s — iv-f iy cos 2sxt—ix{ s ~~ l+2m)sf — iv sin 2s\t 



r\ ml (ve~ 2lsxt ) n L s + n {-ve- 2lsxt ) 



J \nJ (m-r + n)\ L m (-v) 



, (8) 



where L 7 ^ denotes the associated Laguerre polynomial. 
Setting m — 0, we recover Eq. for an initial state 
\a). For an initial state |a, 1), corresponding to minimal 
departure from coherence in this class of PACS, we set 
r = 0, s = 1 and m = 1 in Eq. JSJ to obtain the expec- 
tation value of a(t), and hence the expectation value of 
x(t), denoted by (x(t)) 1 . We get 

(x(t)) 1 = e-^-^^ixoRez^t) +p lm Zl (t)}, (9) 
where 

'2 + ve 2lxt ' 



l + v 



J(2 X t+v sin 2 X t) 



(10) 
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Similarly, for the initial state \a, 2) we find 

(x(t)) 2 = e-^ 1 - cos2 ^[x Rez 2 (t)+p Q lmz 2 (t)}, (11) 
where 



Z2 



( & + Que 2 ^ + v 2 e^ \ (4xt+ „ sin 2%t) 



It is evident that these expressions are already consider- 
ably more involved than that in Eq. Q for to = 0. In 
particular, we note the occurrence, in the exponents, of 
secular terms in t apart from sinusoidal terms involving 
higher harmonics. Similar remarks apply to the expec- 
tation value of p in these cases. The time-dependence of 
the higher moments of x and p are even more involved. 
We do not go into these here. 

In general, the effects of the departure from coherence 
of the initial wave packet on its subsequent dynamics are 
not significant if v m, since the overall exponential 
factor exp [— f(l — cos 2%i)] ensures that the expectation 
values are essentially static except for sudden changes at 
times close to revivals and fractional revivals. This is 
not valid for smaller values of v. In particular, squeezing 
and higher-order squeezing effects depend crucially on 
the precise nature of the initial state, as we shall now 
see. 



III. SQUEEZING AND HIGHER-ORDER 
SQUEEZING 

We begin by summarizing the condition for g th -power 
amplitude-squeezing of a quantum state and outlining 
the squeezing properties exhibited by an initial CS as 
it propagates through a Kerr-like medium [4|. We then 
examine in detail the squeezing properties exhibited by 
an initial PACS as it evolves in the medium, and compare 
its behavior with that of a CS. 

One first defines the two quadrature variables 



Z l — > Z 2 — 



V2 iV2 
The generalized uncertainty principle gives 
(AZ 1 f(AZ 2 ) 2 >i\([Z 1 ,Z 2 ])\ 



(?=1, 2, 3,. 



(13) 



(14) 



where AZi is the standard deviation of Zi , and the ex- 
pectation values refer to those in the state concerned. 
The state is said to be g th -power amplitude-squeezed in 
the variable Z\ if 



(AZ 1 ) 2 < 1 2 \([Z 1 ,Z 2 }] 



(15) 



Amplitude squeezing in Z 2 is similarly defined. We write 
[a q , a^ q ] = F q (N) (this is a certain polynomial of order 
(q— 1) in the number operator N) and define the quantity 



D q {t) 



(AZ^-^{F q (N)) 
k(F g (N)) 



(16) 




D, 



FIG. 1: Plots of D q (^T TCV ) versus v for an initial CS (with 
= 0), for (a) q = 1 and (b) q = 3. (Note the different 
ordinate scales in the two cases.) 



where the time-dependence has been indicated explicitly 
to remind us that the expectation values involved are 
those in the instantaneous state of the system. It is easily 
seen that the state is g th -power amplitude-squeezed in Z\ 
if — 1 < D q < 0. We can rewrite Eq. (|16l) in terms of a q 
and a tlJ as 



D q (t) = 



Re (a 2q ) - 2 (Re (a q ) f + ( ffl t« a q ) 



(F q (N)} 



(17) 



When the initial state is the CS \a), the (time-dependent) 
expectation values on the RHS can be evaluated 4] using 
Eq. ©. 

We are interested, in particular, in examining whether 
fractional revivals are accompanied by any significant de- 
gree of squeezing and higher-order squeezing. For this 
purpose, we focus on D q (t) at the instant t — ir/(2x) = 
iT rcv , corresponding to a 2-sub-packet fractional revival. 
After simplification, wc find that D q vanishes at this in- 
stant for all even values of q, implying that no even-order 
squeezing of the state accompanies this fractional revival. 
For odd values of q, however, we find 



D q (T IOV /2) 



2v q 



(F q (N)) 



sin 



qO), (18) 



where 9 is the argument of a (= v 1 / 2 e 10 ). Thus squeez- 
ing (or higher-order squeezing) occurs at this instant pro- 
vided D q < 0, i.e., |tang0| < e~ 2v . We illustrate this 
in Figs, da) and (b), where Di and D 3 are plotted as 
functions of v for an initial CS with 8 — 0. (We have 
set % = 5 in all the numerical results presented in this 
paper.) We have also plotted D\ and D3 as functions 
of for a fixed value of v (= 0.1) in Figs. Efa) and (b), 
showing how squeezing occurs for certain ranges of the 
argument of a, when D q becomes negative. 

We now show that if the initial state departs even 
marginally from coherence, as in a PACS with a small 
value of m, these results change significantly. Writing 
D q as Dq when the expectation values in Eq. (jT7|l are 
evaluated for an initial state \a, to), we find the following 
general result: setting A = ^L m (—v) (F q (N)}, we have 
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AD^tt) = e -^(i-™s4 xg t) y> ( m + 2 0\ CQS / 2 r 2m + 2 n + q- l) X qt + vsmA X qt - 2q9 

\n + 2q J n! V 

) — j-cos f(g-l + 2m + 2ri)xgt + ^ sin 2xgi — g# ) > 



n + 2q J n\ 

2 e -2^(i-cos2x«3t) fJ^ f I + q\ v n+q 



L m {-v) 



+ y 7 \7 v n L rl (—v), (19) 



where rt m j n = Max(0, 5 — to). This generalizes the ex- 
pression for D q obtained 4] for a coherent state, which 
corresponds to m = 0. 

As before, we examine D ( q m) at t = ?r/(2x) for the 
possibility of squeezing and higher-order squeezing. The 
foregoing expression reduces at this instant of time to 

AD q ^(T rcv /2) = L%(-p) cos 2q9 

cosV 

+ E ^-^T^(-^)- (20) 

^ In/ (to- q + n)l mK ' y ' 

n— n min v 7 v ' 



0.15 




2 1 2 3 4 5 1 2 3 4 5 



V V 

FIG. 3: Plots of D 2 (^T ICV ) versus v (with = 0), for an initial 
PACS with (a) m = 1 and (b) m = 10. 



We use this to analyze various cases numerically. In con- 
trast to what happens for an initial CS, it turns out 
that there is no odd-power amplitude-squeezing for an 
initial PACS. Even-power amplitude-squeezing does oc- 
cur, though, for sufficiently large values of v. This is 
illustrated in Figs. EUa) an d (b), which show the range 
of v for which D2 falls below zero. Figures a) and (b) 
depict the variation of D2 with the phase angle 9 for a 
fixed value of v. 

Turning to the extent of squeezing as a function of 
time, Fig. El compares the temporal variation of the 
standard deviation Ax for an initial CS (to = 0) and 
an initial PACS (to = 1), with a = 1 (hence 9 = 0). 
The horizontal dashed line demarcates the level below 
which the state is squeezed. It is evident that squeezing 
in the rigorous sense accompanies the fractional revival 



at t = ^T rev when the initial state is a CS. This feature 
is suppressed when it is a PACS, although Ax does dip 
down considerably around this fractional revival. 

We have focused on g th -power amplitude-squeezing at 
fractional revivals, as this turns out to provide rather 
more discriminatory signatures of higher-order squeezing 
effects than the other alternative, namely, Hong-Mandel 
squeezing However, a few remarks on the latter are 
in order here. The relevant variables in the case of Hong- 
Mandel squeezing are (a + a^) q )/^/2 and (a — ) q / (i^/2) , 
i.e., essentially the q th powers of x and p. For q = 1, of 
course, Hong-Mandel squeezing is the same as amplitude- 
squeezing, but the two kinds of squeezing differ for q > 2. 
Figs. Eta) and (b) show how ((x— (x)) 4 ), the fourth 
moment of x about its mean value, varies over a revival 
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FIG. 4: Plots of D 2 (±T ICV ) versus 9 (with v = 2), for an initial 
PACS with (a) m = 1 and (b) m = 10. 




0.42 0.46 0.5 0.54 0.58 
t 

FIG. 5: Ax versus time in units of T rev , in the case xo = 
\/2,po = 0. 



period for an initial CS and PACS, respectively. The hor- 
izontal dotted lines indicate the bound on ((x — (a;)) 4 ), 
below which fourth-order Hong-Mandel squeezing occurs 
in this quadrature. An initial CS exhibits such squeezing 
near revivals, and comes close to doing so near the frac- 
tional revival at ^T rcv , but does not actually do so. An 
initial PACS does not display such higher-order squeez- 
ing at any time, although ((x — (x)) 4 ) attains its lowest 
value at revival times. However, fractional revivals are 
marked by rapid oscillations of ((x — (a;)) 4 ), these being 
most pronounced around the 2-sub-packet fractional re- 
vival. These features are enhanced further in the case of 
initial states with larger values of m. 



IV. THE WIGNER FUNCTION AND THE 
NON-CLASSICALITY INDICATOR 

Finally, we examine the Wigner functions correspond- 
ing to the wave packets at instants of fractional revivals, 
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FIG. 6: Plots of ((a; — (x)) 4 ) versus time in units of T rcv for 
initial states (a) \a) and (b) \a, 5), with v = 1. Note the very 
different ordinate scales in the two cases. 



to quantify non-classical behavior during their time evo- 
lution. It is well known that the "extent" to which the 
Wigner function becomes negative (as a function of its 
complex argument) is an indicator of the non-classicality 
of the state concerned. The normalized Wigner function 
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FIG. 7: Plots of the Wigner function corresponding to an 
initial state |q) with a — 1, at (a) t — ^T rcv and (b) t 
Here, and in the succeeding figures, /3j = Re j3 , @2 = Im /3 



W(P; t) (where /? € C) is given by [I 
2 



it 

3 J rev 



W(J3; t) 



Re { E (-1/(2 -^^(^/n!) 1 / 2 



x(2/3) n - i p/„(t)^r i (4|/3| i 



(21) 



where pi n (t) is the density matrix element corresponding 
to the state at time t in the oscillator Fock basis. For an 
initial coherent state la) we have the standard result 



Pln{0) 



a* n a 1 



-M 2 



leading to the well-known expression 



W{0] 0) = -e- 2 \ a -^' 



(22) 



(23) 



for the Wigner function. This is positive definite every- 
where in the complex /3-plane, justifying the appellation 
"classical" for an oscillator coherent state |a). The time 
evolution of the Wigner function under the Hamiltonian 
of Eq. may be computed readily using the represen- 
tation p{t) — J2Tn=o P ln W 10 ( n l f° r t ne density matrix. 
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FIG. 8: Plots of the Wigner function corresponding to an 
initial state \a, 1) with a = 1, at (a) t — 0, (b) t = |T rcv and 
(c) t 
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Figures [7f a) and (b) show how the Wigner function for 
this initial CS behaves at the instants of the 2-sub-packet 
and 3-sub-packet fractional revivals, respectively. 

For the initial photon-added coherent state |a, m), one 
finds 



a l-m a *n-m/J^J 



Pln(0)= , T f , Tl 777 xT- (24) 

ml L m (— v) [I — my [n — my 

Correspondingly, the Wigner function at t = can be 
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expressed in the closed form 

W(J3;0)= T ( > L m {\2(3-a\ 2 )e- 2 \ a -P\ . (25) 

It is to be noted that this is no longer positive definite for 
all complex f3, reflecting the fact that this initial state is 
no longer "totally" classical, as it departs from perfect co- 
herence by the photons that have been "added" to \a) to 



produce the PACS. Figures |Ht a), (b) and (c) are plots of 
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FIG. 10: Plot of 8 versus time in units of T rev . 

W(P ; t) for the case m = 1 at t = 0, \T rcv and |T rov , re- 
spectively. The corresponding plots for the case m = 10 
are shown in Figs. EJa), (b) and (c). With increasing 
m, the oscillations of the Wigner function in the /3-plane 
between positive and negative values become more pro- 
nounced, and the region of non-classicality becomes more 
extensive. 

To get an idea of the degree of non-classicality as a 
continuously varying function of time for each of the dif- 
ferent initial states we have considered, it is instructive 
to consider the non-negative quantity S defined asp^| 

6(t) = Jd 2 (3{\W((3-t)\-W((3;t)^ 

= fd 2 0\W(Jh;t)\-l. (26) 

The larger the value of 6, the greater is the extent of non- 
classicality of the state concerned, although of course 5 
alone does not give a complete picture of the oscillations 
of the Wigner function. In Fig. ^|we have plotted 5 ver- 
sus t for initial states \a) , \a, 1), and \a, 10) where a has 
been set equal to unity. It is clear that in the interval 
between t = and t — T rcv , S is least at the 2-sub- 
packet fractional revival, followed by its values at the 
3-sub-packet and 4-sub-packet fractional revivals. This 
feature becomes increasingly prominent for larger values 
of to, showing that the extent of non-classicality also in- 
creases with m, as expected. 
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